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SUMMARY 

An expression for the diffusion constant for self -diffusion in 
metals is derived, based on the assumption that self -diffusion occurs 
by the vacancy mechanism. The factors determining A in the 

expression D = Ae"^/® 1 , where D is the diffusion constant, Q is 
the energy of activation, R is the gas constant, and T is the 
temperature, are found to be the characteristic frequency of vibra- 
tion of the atoms, the crystal structure, the lattice constant, the 
heat of fusion, the temperature, and the activation energy for an 
atomic jump. A comparison of the calculated values of A with the 
limited experimental data available shows that the calculated values 
are of .the - correct order of magnitude. 


INTRODUCTION 

Phenomena such as age -hardening, annealing, and order-disorder 
transformation depend on the ability of atoms to migrate through an 
atomic lattice. It is therefore of great importance to understand 
the factors that determine the rate at which such migration occurs. 
One method of determining the movements of atoms within a lattice is 
to study the rate at which metals diffuse into each other under 
controlled conditions. A special case of diffusion in metals is the 
case of self -diffusion. 

The amount of material diffusing across a surface because of a 
difference in concentrations of a material is proportional to the 
concentration gradient at the surface, to the area of the surface, 
and to the time during which the process continues. The proportion- 
alities can be converted to equalities by the introduction of a 
factor called the diffusion constant, which has the dimensions area 
per time . This diffusion constant D is not constant for a given 
material but depends on the temperature; the relation is given 
(reference 1) by 
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D = Ae" Q / Rr (1) 

where the coefficient A and the energy of activation Q are 
practically constant for large temperature intervals; R is the gas 
constant and T is the temperature in degrees Kelvin. In order to 
have a better understanding of the diffusion process, it is desirable 
to know the properties of the metals that determine the diffusion 
constant and consequently the rate of diffusion. 

A number of generally known expressions for the diffusion con- 
stants for metals in the solid state exist (reference 2, p. 305). 

The Dushman -Langmuir expression (table I) , which was derived by 
analogy with an earlier expression for the reaction-velocity con- 
stant (reference 2, p. 298), indicates that A depends on the 
lattice constant and on the energy of activation Q; Bradley's 
expression, which was derived on the basis of kinetic theory and on 
the assumption that the distribution of energies among the atoms 
was a Maxwell -Boltzmann distribution/ indicates that A depends on 
the lattice constant, on the frequency of vibration of the atoms, 
on Q, and on the temperature; Cichocki's expression, which was 
derived on a probability basis and on the assumption that the atoms 
become ionized in the process of making a Jump, indicates that A 
depends on the melting temperature a n d on atomic weight ar^i volume . 

More recently, on the basis of rate-process theory, Eyring has 
developed an expression for D (reference 3, p. 538) in which A 
depends on atomic volume, temperature, and the Debye characteristic 
temperature. The preceding expressions were derived for the general 
case of diffusion and as shown in table I only the Dushman -Langmuir 
formula yields results of the correct order of magnitude for the 
special case of self -diffusion. The Dushman -Langmuir formula is, 
however, semiempirical, having been derived without the use of any 
assumptions as to the mechanism by which diffusion occurs and con- 
sequently it does not explicitly show the properties of the metals 
that most directly affect the diffusion constant. 

Although many expressions exist for the diffusion constant, the 
nechanism of diffusion suggested by Frenkel (reference 2, p. 293) is 
gener ally accepted. According to the Frenkel theory, diffusion is 
the result of atoms Jumping into vacant sites, or holes, in the lattice. 
A vacancy is generated in the lattice when an atom Jumps from its 
normal site to an interstitial position, or when an atom escapes 
from its normal site at some internal surface , without occupying 
an interstitial position. Schottky has shown (reference 4, p. 69) 
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that the first process is more likely in a lattice consisting of 
atoms whose radii differ greatly from each other; whereas the second 
process is more likely if the atoms have approximately equal radii. 

In- a theoretical study made at the NACA Lewis laboratory, only 
self -diffusion of metals was considered. 

SYMBOLS 

The following symbols are used in this report: 

A constant 

D diffusion constant 

d average length of atomic jumps 

E f heat of fusion 

E_ heat of sublimation at 0° K 

E Q average energy of vacancy formation 
average thermal energy 

© 

Ei average energy required to make one Jump 

F fraction of atoms having energy equal to or greater than E^ 
J a number of Jumps per second for atom 

J „„ marl nnim possible number of Jumps per second for vacancy 
max 

J 0 number of Jumps per second for vacancy 
m mass of atom 

N number of atoms per unit "cell" 

N a number of atoms per unit volume 

N 0 number of vacancies per unit volume 

Q energy of activation 



4 NACA TN No. 1856 

R gas constant 

T temperature 

t time 

t 0 time required for vacancy to make one Jump 
V velocity of Jumping atom 

v instantaneous velocity of vibrating atom 
Vmax maximum velocity of vibrating atom 
W number of nearest neighbors in lattice 

x displacement of vibrating atom 

xq amplitude of vibrating atom 

v frequency of vibration 

ANALYSIS 

The formation of a vacancy at an internal surface has one of two 
•possible results. The vacancy may either be reoccupied by the atom 
that has Just vacated it, in which case no diffusion occurs, or it 
may be occupied by a different neighboring atom, which leads to dif- 
fusion. It is shown in reference 3 (p. 519) that if J a is the 
average rate at which each atom Jumps, the diffusion constant D is 
given by ‘ 

D = J a d 2 ( 2 ) 

where d is the average length of the Jumps, the Jumps having only 
one degree of freedom. For three degrees of freedom, the equation 
becomes 


D = 



(3) ' 


It will now be shown that J a can be expressed in terms of the num- 
ber of vacancies per unit volume in the lattice • N 0 and the average 
number of Jumps per second made by each vacancy J G . 
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Only atoms adjacent to a vacancy can make jumps . Over a long 
period, however, each atom will have had the same number of oppor- 
tunities for Jumping, so that the average number of jumps per second 
will be the same for all atoms. If a unit volume of the metal under 
consideration is divided into N 0 equal volumes, or 'cells, each 
cell will, on the average, contain one vacancy. If N is the number 
of atoms in a cell, the average number of jumps for any atom in a 
given cell will be l/R times the total number of Jumps. Therefore 
the, numb er of jumps per second for an atom is l/N times the number 
of Jumps per second for a vacancy. In symbol form, 

J a = Jo/n 

The number of atoms N in a cell is given by 

N = Na/N 0 (5) 


where N a is the number of atoms per unit volume . Substitution of 
equation (5) into equation (4) and then of equation (4) into equa- 
tion (3) yields 


a ££o d 2 
3N- . 


( 6 ) 


The value of N c depends on the number of atoms having suffi- 
cient energy E 0 -to escape from their normal sites. For a lattice 
containing vacancies but no interstitial atoms, the total number of 
vacancies per unit volume, when the lattice is in equilibrium, is 
given (reference 4, pp. 63 and 65) by 


N 0 = N a e“ E g/ 2Rr (7) 

where Eg is the energy of sublimation at absolute zero. Because it 
is desirable to know N 0 as a function of E 0 rather than as a 
function of Eg, it is necessary to find a relation between E 0 and 

E S* 

• / 

The following diagram is a schematic representation of an 
internal surface . A vacancy is formed whenever one of the atoms 
acquires sufficient energy to jump from its original position to some 
other equilibrium position. Some atoms, such as atom 1, are so 
located that they can reach a new equilibrium position without appre- 
ciably changing their distance from their nearest neighbors. The 
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energy required for such motion is equal to the heat of fusion. An 
atom such as atom 2, however, must jump a large distance in order to 
break the bonds of its nearest neighbors. It is assumed that the 
distance atom 2 has to jump is large enough that the energy necessary 
for the jump is equal to the sublimation energy. Although a vacancy 
is more likely to form at a place 
such as the site occupied by 
atom 1 because less energy is 
required, the vacancy formed at 
such a point will be partly 
blocked by the atom that has 
vacated the site, thus seriously 
interfering with further diffusion. 

As a first approximation, it will 
therefore be assumed that the 
average energy of vacancy formation 
E 0 is the arithmetic average of 
the minimum and maximum energies of 
vacancy formation. The sublimation 
energy at high temperatures is 
approximately equal to the sublimation energy at absolute zero 
(reference 5), and therefore, E g may be used for the maximum energy 
of vacancy formation. That is, 

*0 - | (Eg + Ef) 

Eg = 2E 0 Ejp (8) 

i 

Substitution of equation (8) into equation (7) yields 

■ 0 = » a e E f/ 2Er e' E o/ OT (9) 

In the determination of J D , the assumption is made that, on 
the average, an atom (providing it has sufficient energy) begins its 
jump as soon as a neighboring atom has exceeded its normal amplitude 
in the course of making a jump* Because the atoms are executing 
simple harmonic motion (in three dimensions), the equation of motion 
in the x -direction is 

x » xq sin 2jtvt 

where x is the displacement, xq the amplitude, v the frequency, 
and t the t ime . Differentiation of the preceding equation yields 



MCA TN No. 1856 


7 


dx 

dt 


v = 2«u Xq cos 2«vt 


When t = 0, v = therefore 


7 max 
x 0 “ 2 «v 


( 10 ) 


When a Maxwell -Boltzmann distribution of energies is assumed, the 
average energy for each degree of freedom is ET (reference 6) . 

The total thermal energy E t is 3ET and therefore the average 
energy for each degree of freedom is l/3 Rfc. The maximum velocity 
v max can then determined as follows: 

\ ^max 2 - !/ 3 Et 


or 

^ Et/3m - (11) 

Substitution of equation. (11) into equation (10) yields 

\]2E t /Zm 
x 0 = Znv 

Let Ei be the energy necessary to make a Jump. Then the 
velocity Y of the atom while it has the energy Ei can be deter- 
mined from the equation 

\ mV 2 = Ei 


or 

Y = hJzE^/m 

If it is assumed that the velocity of the atom remains constant 
while the atom is moving through the distance Xq (xq is only a small 
fraction - about l/25 - of the entire Jump), the time t 0 required 
to traverse the distance xq is 
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t 0 = Xq/v 

^/2E t /3m 

X - 

° ^2Ei/m 

*0 = 2*v 

If t 0 is considered to "be the time required for the vacancy to 
make one Jump, the maximum number of Jumps per second J ma-sr that 
a given vacancy can make is 


Jmax ~ 4- 


2*i> 


'max 


fa /SB] 


( 12 ) 


But a Jump occurs only when one of the atoms near a vacancy has an 
amount of energy equal to or greater than E]_. For a Maxwell - 
Boltzmann distribution of energies, the fraction of the atoms having 
energy equal to or greater than E^ is 


F = 


-Et/ET 

e 


(13) 


-E /ET 

Because at all times el' of the atoms have sufficient energy 
to make a Jump, it is also true that each atom has sufficient energy 
-El ART ^ 

e x/ of the time. If each vacancy had one nearest neighbor, 
the number of Jumps per second would be given by FJ nwv . But if 
each vacancy has, for instance, W nearest neighbors, the number of 
Jumps per second J D is given very nearly by 


Jo = ^max ( 14 ) 

When the values of F and Jmax as given by equations (13) and (12) 
respectively, are substituted into equation (14), there results 


J 0 = 2rtV W fiE 1 /E t e ~ El ^ 


( 15 ) 


Substitution of equations (15) and (9) 


into equation (6) yields 
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D = §SL W * 2 e Ef/2BT e" (E ° + El)/Er (16) 

Consequently, from a comparison of equations (l) and (16) 

A = W */3 Eq/Et d 2 & E f/ 2ET (17) 


COMPARISON OF CALCULATED AND EXPERIMENTAL DATA 

All the data available for metals with symmetric crystal struc- 
ture are for face-centered cubic lattices. The number of nearest 
neighbors is the same for all materials of this type and is equal 
to 12. When all the constants of equation (17) are combined, A 
becomes 

A = 43 . 6Vd 2 ^Ei/Et e E f/ 2Rr (18) 

The average length of the Jumps is assumed to be equal to the lattice 
constant. Equation (18) can be applied to determine the value of A 

if V E i/ E t can be approximated, inasmuch as all the other parameters 
are known. 

On the basis of existing data (reference 2, p. 275) the maximum 
value of Q is about 60,000 calories per mole, that is, 

E 0 + Eq S 60,000 

Theoretically the s mal lest value that E 0 can have is the heat of 
fusion Ef, which for copper, gold, and silver is approximately 
3000 calories per mole. The smallest value for E 0 yields the 
largest theoretically possible value for Eq . That is, if E 0 = Ef 

Ef + Eq s 60,000 

Eq S 60,000 - 3000 S 57,000 

fife 2^57,000/Et 

and 


E t * 3RP S 6T 
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The center of the temperature range in which experiments are conducted 
is approximately 1000° K (for copper, gold, and silver) . At this 
temperature, the ma-rlmmm theoretically possible value for 

^E 1 /E t = -s/57, 000/6000 = 3.08. 

It is more reasonable to assume that 

E o = E 1 

For this case, 

ZE ± = 60,000 
E 1 = 30,000 

( 

and 

^Ex/E-fc = \J 30,000/ 6000 S 2.24 

For lead, the center of the temperature range is about 500° K and Q 
is about 28,000 calories per mole. Therefore, 

E 2 £ 14,000 


and 

^E 1 /E t = /^/l4,000/3000 £ 2 .16 

From the preceding considerations it can be seen that a/e^/e^ is to 
a large extent independent of material as well as o f temp erature, 
having a value approximately between 2 and 3. If ^E^/E-t is set 
equal to 2.5, equation (18) becomes 

A = 109 vd 2 e Ef / 2Rr (19) 

Although equation (19) shows that A is a function of the temperature , 
the calculated value of A is only about one and one-half times as 
great at 750° K as it is at 1250° K. Inasmuch as the experimental 
values of A, which are independent of temperature, show a greater 
variation, A will only be ev al uated at 1000® K for gold, coppex* , and 
silver, and at 500° K for lead. 

The value of A can be calculated for the metals for which self- 
diffusion data exist by using the data in table II • The experimental 
and calculated values of A are contained in table I . 
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Another comparison can he made with previous calculations. On 
the basis of the random-walk theory, Weinbaum (reference 7) has cal- 
culated that the average number of jumps per second J a for the 
atoms of a silver lattice at a temperature of 1149° K (876° C) is 

6.3 x 10®. Equations (4) and (5) can also be used to calculate the 
number of jumps per second under the aforementioned conditions. 

a N a 

Substituting equations (15) and (9) for J Q and N 0 , respectively, 
and combining constants yields 

J a = 324 v e E f/ 2OT e- (E ° + E 1 )/ET 

For silver, Q = 45,950 calories per mole (reference 8). At 
1149° K 


J a = 5.8 X 10 6 
SUMMARY OF ANALYSIS 

An expression for the diffusion constant for self -diffusion 
has been derived and has been found to give results that are of the 
right order of magnitude. The coefficient A in the expression 

D = Ae - ^^, where D is the diffusion constant, Q is the energy 
of activation, R is the gas constant, and T is the temperature, 
is determined by the characteristic frequency of vibrations of the 
atoms, the crystal structure, the lattice constant, the heat of 
fusion, the temperature, and the activation energy for an atomic 
Jump. More data are needed, however, before definite conclusions 
can be drawn regarding the applicability of the expression. 


Lewis Flight Propulsion Laboratory, ' 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, February 2, 1949. 
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TABLE I - COMPARISON OP 
VALUES OP 

D = 


CALCULATED AND EXPERIMENTAL 


A FOR EQUATION 
Ao-^ T 


' Metal 

Experimental 
values of A 
( cm /sec ) 

Calculated values of A 
(cm^/sec) 



Reference 

(a) 

(b) 

(c ) 

(d) 

(e) 

(f) 

Lead 

5.1 

6.66 

2 , ?. 274 
9 

0.88 

0.71 

2.0 x io' 2 

5.7 x io* 5 

4.2 x io" 5 

4.5 x 10 5 

Silver 

1.0 

.895 . 

10 

8 

1.6 

0.7k 

.80 

2.6 x io’ 2 
2.8 

6.4 x io' 5 

2.1 x io" 5 

7.8 X 10 ^ 

Copper 

ll 

47 

2 , p. 274 
9 

2.0 

0.78 

.64 

4.2 x 10‘ 2 

4.4 

7.3 x io* 3 

2.5 x xo" 5 


Gold 

2 

126 

410 

9 

2, p. 274 
9 

1.5 

0 . 88 
.88 

1. b 

2.6 x 10* 2 

2.6 

3-2 

5.5 x io* 3 

1.7 x io' 3 

1.2 x IQ 6 



a A =-2252 V^Ej/Et <* 2 e E ^ 2RT (equation ( 17 )). 


A = "Nh d2 ( Dushman-Langmuir ) , where N ia Avogadro's number and 
h is Planck's constant* 

A = T rt~ pd2 (Bradle y>* 

d A = §£ ( van Llempt ) . 

5rr 


e a - 


( V\ kT -39/lj.T 

¥) IT (1 " ° } (Eyrln « ) ' wh8re V is the gram atomic volume, 

k Is the Boltzmann constant, and 9 Is the Debye characteristic temperature, 

f A = 2.45 x lc)i| y /-jjf (Cichocki), where M, V, and T s are atomic weight, 
and mi T s y°see Reference ^11? ten, P erature ’ respectively. For values of M, V. 
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TABLE II - PARAMETERS FOR EQUATION (19) 


Metal 

Ef 

(cal/mole) 

(a) 

V 

(vibrations/ sec ) 

d 

(cm) 

(c) 

Copper 

3110 

b 6.6'X lCP-2 

3.61 X 10-8 

, Gold 

3030 

*3.8 

4.07 

Silver 

2700 

*4.5 

4.08 

Lead 

1220 

*1.8 

4*94 


a Data from reference 5>. p. ,259. 

*Data from reference 12. 
c Data from reference 13. 

^Computed from v = k9/h where k is Boltz- 
mann constant, h is Planck's constant, and 
9, the Debye temperature for gold, is 180 
(reference 5, p. 237). 



